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[Document Name] SPECIFICATION 
[Title of the Invention] Method of ML-EM image reconstruction 

and medical imaging apparatus 

[Claims] 

1. A method of ML-EM image reconstruction for use in connection with a diagnostic 
imaging apparatus that generates projection data comprising: 

(a) cbHectmg projection data including measured emission projection data; 

(b) assuming in initial emission map; 

(c) assuming an initial attenuation map; 

(d) iteratrvely updating the emission map and the atte nu ation map; 

(e) with each iteration, recalculating the emission map by taking a previous 
emission map and adjusting it based upon: CO the measured emission projection data; (n) 
areprojecdon of the previous emission map which is earned out with amuld-dimensional 
projection model; and (m) a reprojectiDn of the attenuation map; and 

(f) with each iteration, recalculating the attenuation map by taking a previous 
attenuation map and adjusting it based upon: (L) the measured emission projection data; and, 
(H). the reprojection of the previous emiiaion map which is carried out with the multi- 
dimensional projecti on model. 

2. A method as claimed in claim 1, wherein the step of collecting projection data 
further includes collecting measured transmission projection data. 

3. A method as claimed in claim 2, wherein the adjustment to the previous attenuation 
map in the recalculating step is also made based upon the measured transmission projection 
data. 
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4. A method as claimed xo claim 2 or claim 3, wherein the step of collecting projection 
data further includes collecting a blank set of transmission projection data which ia 
collected without a subject present in the diagnostic imaging appantus, and the adjustment 
to the previous attenuation map in the recalculating step is also made based upon the blank 
set of transmission projection data. 

5. A method as claimed in any one of claims 2 to 4, wherein the measured transmission 
projection data is truncated. 

6. Medical imaging apparatus for reconstructing image representations of a subject 
being examined therewith comprising: an emission memory (113) which stores measured 
emission projection data from the subject; an image memory (140) which stores emission 
maps; an attenuation map memory (130) which stores attenuation maps; a first projector 
(ISO) which generates emission map projections via forward projection of the emission 
maps from the image memory (140), said projector (150) using a first multi-dimensional 
projection model; a second projector (152) which generates attenuation map projections via 
forward projection of the attenuation maps from the attenuation map memory (130), uiA 
projector (130) using a second multi-dimensional projection model; a first data processor 
(160) which samples the emission memory (112), and the first projector (150), and in 
accordance therewith updates the attenuation maps stored in the attenuation map memory 
(130); a second data processor (170) which samples the emission memory (112), the first 
projector (150), and the second projector (152), and in accordance therewith updates the 
emission maps stored in the image memory (140); and a human-viewable display (310) 
which renders reconstructed images representations of the subject from the emission maps. 
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7. Medical imaging apparatus as daimed in elaim 6, further comprising! ^ transmission 
memory (114) which stores measured transmission projection data from the subject, 
wherein the medical diagnostic imaging apparatus (10) collects said transmission projection 
data. 

8. Medical imaging apparatus as claimed in claim 7, wherein the first data processor 
(160) also samples the transmission memory (114) when it updates the attenuation maps 
stored in the attenuation map memory (130), such that the update is also made in 
accordance therewith. 

9. Medical imaging apparatus as claimed in claim 7 or claim 8, fijithcr comprising: a 
reference memory (120) which stores transmission projection data from a blank run 
performed without the subject present in the medical diagnostic imaging apparatus (10), 
wherein the first data processor (160) also samples the reference memory (120) when it 
updates the attenuation maps stored in the attenuation map memory (130), such that the 
update is also made in accordance therewith, 

10. Medical imaging apparatus as claimed in any one of claims 7 to 9, wherein the 
measured transmission projection data from the subjected stored in the transmission 
memory (114) is truncated. 
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[Detailed Description of the Invention] 
[0001] 

[Technical Field to which the Invention pertains] 
The present invention relates to imaging using image reconstruction of amission data. It 
finds particular application in conjunction with single photon emission computed 
tomography (SPEC!) and nuclear cameras, and will be described with particular reference 
thereto. However, it is to be appreciated that the present mvention is also amenable to other 
like applications including, but not limited to, positron emission tomography (PET) and 

ultrasound tomography. 

[0002] 
[Prior Art] 

Attenuation correction is now an integral part of nuclear medicine image reconstruction. 
Several techniques have been introduced over die years, ranging from simple analytical 
filters to elaborate transmission devices, to obtain an object-dependent attenuation map. 
Often, in analytical correction, assumptions are made that are too restrictive to be realistic. 
Object-dependent attenuation correction based on emission data only is an attractive and 
theoretically sound alternative, but not ready for large scale utilization. Finally, techniques 
based on an a tt e nu ati o n map constructed from measured transmission projection data, while 
a popular approach, suffer from very poor statistical quality. Moreover, current techniques 
for reconstructing the emission data are sensitive to this noise. Additionally, the collected 
transmission projection data is often truncated. 
[0003] 

The approaches presented herein are a departure from previously available techniques. In 
order to appreciate the difference, it is helpful to examine previous techniques. 
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It is usually recognized that iterative methods lend themselves easily to the introduction of 
additional correction ftcton such as attenuation. With reference to FIGURE 1, the buic 
concept behind iterative methods, such as maximum likelihood expectation m aximizat ion 
(ML-EM) reconstruction, is to find the object (Lc, image or emission map) for winch & 
mathematical projection produces results comparable to the set of measured emission 
projections. The ML-EM reconstruction algorithm is one such engine that allows the image 
to be found in an efficient manner. From an initial guess or assumption, a set of projections 
is artificially created by a projector employing a projection model. These projections are 
then compared to the "rear or measured set When certain conditions are met (U, when 
the projections of the emission map are sufficiently dose to the measured emission 
projections), the iterative process stops and the current image is the best possible 
representation of the object, otherwise, the initial guess is -updated and a sew set of 
projections produced. Clearly, the projection model employed is an important part of the 
projection operation. However, an accurate projection model can be complex and detailed. 
[0004] 

With reference to FIGURE 2, one way to improve the "realism" of the projection model or 
the update operations is to include a priori information, ie., what is already known about 
the object (e.g^ its contour, texture, and so forth). In general, this is a powerful tool that 
offers dramatic improvement of the reconstruction. One caveat is, however, that this 
information needs to be teal In other words, if a priori information is true or accurate, it 
helps the reconstruction. On the other hand, inaccuracy tends to inappropriately bias the 
reconstruction or otherwise introduce unwanted arti&cts. Moreover, in practice, it is 
difficult to find non-trivial characteristics of the object to be imaged that ire always 
accurate. 
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[000&] 

With reference to FIGURE 3 , the attenuation map is certainly one element that can help the 
reconstruction. Information pertaining to a t ten u a t i o n characteristic* can be generated 
artificially for simple situations (e.g. uniform attenuation, symmetrical attenuation, etc). 
In general, however, inaccurarica can be introduced by mggesting attenuation feature* that 
are sot tine and/or not object specific 
[0006] 

With reference to FIGURE 4, it is significantly more advantageous when the attenuation 
map is constructed for each object and used directly in the reconstruction. Often, the 
attenuation map is derived from measured transmission projection data. This is the basis, 
of most of the modem non-unifbrm attenuation correction devices. However) the problem 
with this approach is statistics. Typically, in order to define a useful and usable attenuation 
map, many counts are needed. In practice, however, the definition of the attenuation map 
is count-limited, and the inherent noise associated with it is transported into the 
reconstruction of the emission map. Another problem is that often the transmission 
projection data used to derive the attenuation map is truncated thereby resulting in a 
degraded attenuation map. In other words, a poor attenuation map can actually degrade the 
ltmgi» it aims at i mp ro v ing because of a lack of counts and/or a truncated transmission data 
set 

[0007] 

These limitations in the acquisition of the measured attenuation maps have provided the 
motivation for investigation into different techniques. It is generally known that the 
emission data contains information from which theoretically the attenuation map could be 
reconstructed. Such reconstruction methods are referred to as sourceless. In fact, each 
point in the object can be considered as a transmission source, and the observed Intensity 
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it any given point on the detector can be compared with the expected intensity without 
at te nuati on. Unfortunately, due to limitations, such u the Kmhed number of counts in the 
emission map, the relationship between the emission and attenuation maps can only 
partially be established. Moreover, there is an intrinsic inability of the process to 
differentiate between a low activity, low attenuation condition and a Ugh activity, high 
attenuation condition. 

[0008] 

Yet another drawback of previous techniques is that they tend to employ only a linear 

projection model That is to say, they merely account for those events that lie along a 

singular path or fine from the detection bin to the activity. Such projection models are not 

adapted to account fa- nr collimator ^^frrtifln whir h may rites* a particular detector 

bin to detect photons from off-line sources. 
[0009] 

[Problems that the Invention is to solve] 

The present invention contemplates a new technique for providing attenuation correction 

in emission computed tomography. 
[0010] 

[Means for solving the Problems] 

In accordance with one aspect of the present invention, a method of ML-EM image 
reconstruction is provided for use in connection with a diagnostic imaging apparatus that 
generates projection data. The method includes collecting projection data including 
measured emission projection data. Initial emission and attenuation mapa are ass umed a n d 
iteratively updated. With each iteration, the emission map is recalculated by taking a 
previous emission map and adjusting it based upon: © the measured emission projection 
data; 00 * reprqjection of the previous emission map which is carried out with a xmihi- 
cfimensiooal projection model; and, (Si) a reprojection of the attenuation map. Also, with 
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each iteration, the atten uation trap is recalculated hy talcing a pngygftiM ^ttnn?fltion !tt*p tlld 
adjusting it based upon; © the measured emission projection data; and, (H) thereprojeetion 
of the previous emission map which is carded out with the muW-dimensional projection 
m od el. 

[0011] 

In accordance with another aspect of the present invention, a medical diagnostic imaging 

apparatus for reconstructing image representations of a subject bang examined therewith 

is provided. It includes an emission memory which stores measured emission projection 

data from the subject An image memory stores emission maps, and an attenuation map 

memory stores attenuation maps. A first projector generates emission map projections via 

forward projection of the emission maps from the image memory, likewise, a second 

projector generates attenuation map projections via forward projection of the attenuation 

maps from the attenuation map memory. The first and second projectors use first and 

second muM-dimensional projection models, respectively, A first data processor samples 

the emission memory, and the first projector, and in accordance therewith updates the 

attenuation maps stored in the attenuation map memory. A second data processor samples 

the emission memory, the first projector, and the second projector, and in accordance 

therewith updates the emission maps stored in the image memory. Finally, a human* 

viewable display renders reconstructed image representations of the subject from the 

emission maps. 
[0012] 

[Embodiments] 

Ways of canyixjg out the invention will now be described in detail, by way of example, with 
reference to the accompanying drawings, 
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With reference to FIGURE 5, the basic flow of a sourceless ML-EM reconstruction in 
accordance with the present invention is illustrated. A medical diagnostic imaging 
apparatus 10 (described in greater detail below) generates emission projection data. In this 
embodiment, the transmission map is reconstructed from emission data only. Lack of a 
transmission scan reduces the patient's radiation exposure. Moreover, in a preferred 
embodiment a multi-dimensional projection model provides a more accurate ^projection 
of the emission map. With reference to FIGURE 6, the basic flow of a source-assisted 
reconstruction is illustrated which is a further adaptation of the sourceless reconstruction. 
In this embodiment, the medical diagnostic imaging apparatus 10 also generates 
transmission projections. Source-assisted attenuation correction combines the benefit of 
the real transmission measurements to the stability of the analytical sourceless approach. 
In this approach, the measured transmission projections, however statistically poor and/or 
truncated, are used to stabilize the analytical approach and force the attenuation map and 
the emission map to be independent The transmission information and emi ssion counts 
are combined to form an image of the attenuation map. The result is improved quality in 
the attenuation map, because the emission data is used to "assist" in its reconstruction. 
This, in turn, improves the reconstruction of the emission map in terms of quantitative 
accuracy. 
[0013) 

The medical diagnostic imaging apparatus 10 is generally any nuclear medicine scanner that 

* 

produces scintigraphic images, for example, SPECT scanners or PET scanners. An 
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appropriate device is one that detects and records the spatial, temporal, and/or other 
characteristics of emitted photons. 
[0014] 

More specifically, with reference to FIGURE 7, In an exemplary embodiment, the 
£agnostic nuclear Imaging apparatus or scanner 1 0 is a gamma camera or SPECT scanner 
including a subject support 12, such as atahle or couch, which supports and positions a 
subject being examined and/or imaged, such as a phantom or patient, within a subject- 
receiving aperture 18. The subject is injected with one or more radiophannaceutical or 
radioisotopes such that emission radiation is emitted therefrom, A first gantry 14 holds a 
rotating gantry 16 mounted thereto. The rotating gantry 16 defines the subject-receiving 
aperture IS. One or more detector heads 20 are adjustably mounted to the rotating gantry 
16 with varying degrees of freedom of movement Being mounted to the rotating gantry 16, 
the detector heads 20 rotate about the subject-receiving aperture IS (and the subject when 
located therein) along with the rotation ofthe rotating gantry 16. In operation, the detector 
heads 20 are rotated or induced around the subject to monitor radiation from a plurality of 
directions to obtain a plurality of different angular views. 
[0016] , 

Each of the detector heads 20 has a radiation-receiving fkee facing the subject-receiving 
aperture IS that includes a sdntfllition crystal, such as a large doped sodium iodide crystal, 
that emits a flash of light or photons in response to incident radiation. An array of 
photomultiplier tubes receives the fight and converts it into electrical signals. A resolver 
circuit resolves the x, y-coordinates of each flash of light and the energy of the in ci dent 
radiation. That is to say, radiation strikes the scintillation crystal causing the scintillation 
crystal to scintillate, ie., emit Rght photons in response to the radiation. The photons are 
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directed toward the pbotamiiltipfier tubes. Relative outputx of the fhnt/wiyiftjpHf tubes are 
processed and corrected to generate an output signal indicative of (1) a position coordinate 
on the detector head at which each ndiation event is received, and (2) an energy of each 
event Tie energy is used to differentiate between various types of radiation such as 
multiple emission radiation sources, stray and secondary emission radiation, transmission 
nufiatiou, and to eliminate noise. An image representation is defined by the isolation data 
received it each coordinate. The radiation data is then reconstructed into an image 
repres ent a t ion of the region of interest 
[0016] 

Optionally, the detector heads 20 in dude mechanical coillimators 24 removably mounted 
on the radiation receiving fkcei of the detector heads 20. The collimators 24 preferably 
include an any or grid of lead or otherwise radiation-absorbent vanea which restrict the 
detector heads 20 from receiving radiation not travelling along selected rays in accordance 
with the data type being collected (te., parallel beam, fin beam, and/or cone beam). 
Alternately, for example, when operating in a coincidence mode, the collimators 24 may 
be omitted. 
[0017] 

One or more radiation sources 30 are mounted across the subject-receiving aperture 18 from 
the detector heads 20. Optionally, they are mounted between the detector heads 20 or to the 
radiation receiving faces of opposing detector heads 20 such that transmission radiation 
from the radiation sources 30 is directed toward and received by corresponding detector 
heads 20 on an opposite aide of the subject-receiving aperture IS. In a preferred 
embodiment, the collimators 24 employed on the detector heads 20, in effect, collimate the 
transmission radiation. Hat is to say, the collimators 24 restrict the detector heads 20 from 
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receiving those portions of transmission radiation not tntvelling along rays normal (pat 
parallel beam configurations) to the radiation receiving frees of the detector heads 20. 
Alternately, other animation geometries are employed and/or the ccffimatiou may take 
place at the source. 
[0018] 

In a preferred embodiment, the radiation sources 30 are moving point sources each 
traversing the length of the respective detector heads 20 to which they correspond 
Alternately, line sources are employed such as, eg., thin steel tubes filled with 
radionuclides and sealed at their ends. In other embodiments, the radiation source 30 are 
bar sources, point sources, flat rectangular sources, disk sources, flood sources, a tube or 
vessel filled with radionuclides, or active radiation generators such as x-ray tubes. 
Alternately, one or more point sources of transmission radiation may be utilized. 
[0019] 

With reference again to FIGURE 7, the running of an imaging operation includes an 
iterative ML-EM reconstruction technique wherein emission projection data is 
reconstructed via an image processor 100 into an image representation of the distribution 
of ratEoactive materia] in the patient Of course, the reconstruction technique implemented 
varies slightly in accordance with the types of radiation collected, the types of collimators 
used (i.e. fan, cone, parallel beam, and/or other modes), and the diagnostic imaging 
apparatus employed. In any case, emission radiation from the patient is received, for 
•ample, by the detector beads 20 of the SPBCT scanner, and measured emission projection 
data is generated. The measured emission projection data normally contains maccuracics 
caused by varying absorption or attenuation characteristics of the patient's anatomy. 
Optionally, a transmission scan is also performed such that transmission radiation from one 



(25) 2000-180550 

or more of the transmission radiation sources 30 is also received, for example by the 
detector heads 20 of the SPECT scanner, and measured transmission projection data is 
generated. However, the measured transmission projection data is often noisy and/or 
truncated. Where a transmission scan is performed, a sorter 200 sorts the measured 
emission and transmission projection data on the basis of their relative energies. The data 
is stored in a measured projection view memory U0, more specifically, in a corresponding 
emission memory 112 and transmission memory 114. 
100201 

In those instances where transmission projection data is collected, a pre-acan blank run is 
performed to collect a baseline or reference with which to evaluate and/or compare the 
measured transmission projection data against. In the blank run, the scanner is used to 
perform a transmission scan without the subject presort therein such that a blank set of 
transmission projection data is collected. In a preferred embodiment, this data is stored in 
a reference memory 120. 
[0021] 

In preparation for the first iteration of the reconstruction process, an attenuation map 
memory 130 and an image memory 140 are initialized by loading them with assumed or 
first estimates of the a tte n uation map and emission map respectively, The first estimate for 
the attenuation map is optionally characterized by a uniform attenuation value inside a 
predetermined contour which contains the subject and zero outside the contour. Likewise, 
the first estimate for the emission map is optionally characterized by a uniform value inside 
the contour and zero outside. Alternately, the availability of additional a priori information 
allows for more accurate first estimates. 
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[0022] 

With each iteration of the reconstruction process, a first projector 150 reprojects and/or 
forward projects the estimated emissioii nap stored in the image memory 140 to obtain 
projection views thereof Likewise, a second projector 152 reprojects and/or forward 
projects the attenuation map stored in the attenuation map memory 13 0 to obtain projection 
views thereof The projections ire performed using multi-dmeniioDal projection models. 
That is to say, the multi-dimensional projection models account for radiation from multiple 
directions, as opposed to merely linear models that only account for radiation travelling 
along a singular direction. In a preferred embodiment, the projection models account for 
any one or more of the following: scatter, the geometric point response (resolution) 
sssodated with the imaging apparatus, randomness corrections which account for random 
simultaneous photon detection not resulting from the same annihilation event in PET 
applications, and/or cross-talk due to changing photon energies. In at least one 
embodiment, the projection models used by the first and second projectors ISO and 152 are 
not the same. Alternately, for some applications, identical projection models are employed. 

A first data processor 160 updates foe attenuation maps in the attenuation map memory 130 
with each iteration, and a second date processor 170 updates the emission maps in the 
image memory 140 with each iteration. 
[0023] 

The first data processor 160 recalculates the attenuation map by taking a previously stored 
attenuation map from the attenuation map memory 130 and adjusting it based upon a 
number of factors. In a preferred embodiment, employing so called sourceless 
reconstruction, the factors include: (3) the measured emission projection data from the 
emission memory 112; and, 09 a reprojection of the previous emission map stored in the 
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linage memory 140 which U xeprpjected via the projector 150. More specially, the first 
data processor 160 samples the emission memory 112 and the output projection view* flora 
the projector 150 to complete each update of the attenuation map which is then reloaded 
into the attenuation map memory 130. For sourceiess reconstruction, measured 
transmission projection data is not utilized. As such, it is preferred that transmission 
projection data not be collected to thereby minimize radiation exposure and reduce 
mechanical complexity. In an alternate embodiment, employing so called source-assisted 
reconstruction, in addition to the above, measured transmission projection data (optionally 
truncated) is also incorporated with each update of the attenuation map. In this 
embodiment, the measured transmission projections are also sampled from tbe transmission 
memory 114 by the first data processor 160 for eacA update or recalculation. In addition, 
baseline or reference transmission projections from tbe reference memory 120 are also 
taken into account A detailed account of tbe first data processor's operation is given below 
in more mathematical terms. 
[0024] 

In both sourceiess and source-assisted reconstruction, the second data processor 170 
updates the emiss ion map based on the same fcctors. These fictors include:' (i) the 
measured emission projection data; 00 a reprojection of the previous emission map; sad, 
Oil) a itprojection of the attenuation map. More specifically, with each iteration, the second 
data processor 170 samples the emission memory 112, the first projector 150, and tbe 
second projector 152. Based upon the data collected, the second data processor 170 
recalculates the emission map and reloads the updated emission map back into tbe image 
memory 140. Again, a detailed account of the second data processor's operation is given 
below m more mathematical terms. 
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[0026] 

Upon completion of the final iteration, a video processor 300 extracts or otherwise receive, 
data from the image memory Uo which represents a reconstruction of the radionuclide 
di5tribution within the region of interest The video processor 300 then formats it fato 
image representations of the fubject for viewing on * huaan-viewhle display, such as a 
video monitor 310 or other appropriate rendering device. 
[0026] 

In more mathematical terms, the operation of the image processor 100 is as fellow,. Tie 
fenrfval^EMreconrtn^ Ith „ 
arobu^inreoonttruct^ 

baaed on the aohition of h«ar.y«ems. ^important feature of the MI^EM algorithm is 
that given certain mM restrictions appropriate for nuclear medicine emission imaging, ft 
ia applicable to any linear system of aquations, whether it be over- or undernietcammed, 
consistent or ir«msisteDt. In any event, after a number of iteration, the ML-EM algorithm 
converge, toasolutioi: that is me mostly 
[0027] 

One of the appealing qualm* of this iterative algorithm is that a general probability model 
ofdetection iiuaed. The underlying radionuclide distribution is continuous, but the signal 
or projection is detected at a discrete number of bins, or sampling location,. In what 
follows, P, represents the total detected photons In the i* detector location (La., the 
measured emission projection data). There are M detector locations. Also, the radionuclide 
distribution is represented as a discrete collection of basis function,. The basis function for 
the enuarion dtaribution is the vo»L F k represent the average of the continuous 
wdioiwcEoe distribute There are K image voxels. The transfer 

matrix codfideat (Le., the projection modeO is the probability that a photon emitted in 



(29) 



2000-180550 



voxd k per unit time is detected in detector location i. Aj each measurement can be noisy, 
the noise associated with each measurement is denoted by n,. Generally however, enough 
is known about the imaging system to be able to accurately model the Jtatirical properties 
of the noise. Hence, the total photon count at the i* detector location is the sum of all 
voxels contributing to that , detector location. This is mathematically expressed as: 



where At| is the imaging time for the i* projection. 
[0028] 

This linear system of equations is then the basic imaging modeL Denoting £™ (Le^ the 
k* voxel of the emission map) as the estimate of the activity F k in voxel k on the h* 
iteration oftheMLUM reconstruction, the M^EM estimate for the (h+l>^ 

br 
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With this method scatter, resolution, snd sttmu stion can be incorporated into the model ty 
However, the inclusion of attenuation and scatter in the transfer matrix employs a priori 
knowledge of the attenuation map. 
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[0029] 

Sonrcelesi Reconstruction for PET Embodiment 

The lade of an attenuation map for emission imaging in this case leads to a new formulation 
that now includes the unknown attenuation map. The forward problem ia now stated as: 

^i-At^a^^-^g Cil u^to J# . \ < 3 ), 

where ia the linear attenuation coefficient of the j* pixel of the attenuation map, and Cj 
ia the effective attenuation length of fij. The set of pixels N M is given by: N M ** {pixel j | 
% + 0, for aouiccnletector pair (F b PJ}. Geometrically, the aet N^is the aet of pbcels 
intersecting and covering the line of response i that contains the source pixel k. The 
estimate (f k , iij) of th e unknowns (F fc , ty) is given fay: 

for the emission map which is updated via the second data processor 170; and 

■T" - «/" - — ^ v . ' <s) , • 

for the attenuation map which is updated via the first data processor 160. That is to say, in 
this embodiment (le., soundest reconstruction with PET data), the second data processor 
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170 tad the first data processor 160 carry out the fanner tad later computations from 
equations (4) and (5), respectively. 

[0030] 

Source! Reconstruction for SPECT Embodiment 

For SPECT, a similar set of algorithm! ii developed. The forward problem hi this cue is 
again: 



However, for SPECT, the wc N M a defined differently as illustrated in FIGURE 8, The set 
N M is the sec of pixels that cover the attenuation line integral from the source pixel k to the 
projection bin i. This leads to a variation of the algorithm. In this embodiment, the 
emission update is given by: 



and, the transmission image update is given by: 
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As before, in this embodimeot (L e^ sourceless reconstruction with SPECT data, the second 
data processor 170 and the first data processor 160 cany out the former and later 
computations from equations (7) and (8), respectively. 
[0031] 

Source-Assisted Reconstruction for PET Embodiment 

Source-assisted reconstruction utilizes the measured transmission data, represented by T v 
and blank scan, represented by E^ to stabilize the sourceless reconstruction. It 
accomplishes this by incorporating the transmission scan as a priori information in the 
likelihood probability Amotion. The result Is a different update for the transmission map 
only. For PET, the transmission update is now given by; 



(9), 



in which Atj is the scan time per projection for flic emission data, and is a normalization 
constant that reflects the ratio of the emission scan-time to the transmission scan time. The 
coefficients 6 jn measure the contribution erf the j* voxel of the attenuation map to the n* 
transmission measurement and are not always the same as the coefficients c^ defined above. 
They depend on the particular embodiment of the scanner employed. 
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[0032] 

Source-Ait isted Reconstruction for SPECT Embodiment 

The update rule for the emission map in source-assisted SPECT is the same as the update 
rule for lourcdcis SPECT. Only the transmission nup update rule is changed to the 
following: 



Optionally, the operations of the first and second data processors 160 and 170, ts well as 

the other components, are implemented as hardware^ software, or combinations of hardwire 

and software configurations. 
[0033] 

[Advantages of the Invention] 

One advantage of the embodiments described above of source-assisted attenuation' 
correction for emission computed tomography! is improved attenuation maps for ML-EM 
reconstruction. Another advantage is reduced patient exposure to radiation via sourceleas 
ML-EM reconstruction. Yet another advantage is multidimensional projection modelling 
which more accurately approximates actual conditions, Another advantage is that there is 
less likelihood of noise being introduced into the final reconstruction via the attenuation 
map. A further advantage is that accurate reconstruction is possible even with truncated 
transmission data. 
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[0034] 

[Brief Description of the Drawings] 

[FIG, 1] lit flow diagram illustrating the operation of a typical MUEM reconstruction 
method in accordance with prior art techniques; 

[FIG. 2] is a flow diagram illustrating the operation of a typical ML-EM reconstructs 
method employing a priori information in accordance with prior art techniques; 

[FIG. 3] is a flow diagram illustrating the operation of a typical ML-EM reconstruction 
method employing an artificially generated attenuation map in accordance with prior ait 
techniques; 

[FIG. 4] is a flow diagram illustrating the operation of a typical ML-EM reconstruction 
method employing an attenu a ti on map derived from measured transmission projection data 
in accordance with prior art techniques; 

[FIG. 6] is a flow diagram illustrating aourcelesa ML-EM reconstruction in accordance 
with aspects of the present invention; 

[FIG. 6]; la a flow diagram illustrating source-assisted ML-EM reconstruction in 
accordance with aspects of the present invention; 

[FIG. 7] is a diagrammatic illustration of a diagnostic nuclear imaging apparatus 
employing an image processor which carries out ML-EM reconstruction is accordance with 
aspects of the present invention; and 

[FIG. 8] is a graphical illustration showing the set of pixels employed in the 
reconstruction of SFBCT data in accordance with aspects of the present invention. 
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FIG. 8 
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[Document Name] ABSTRACT 

[Abstract] 

[Problems] 

There is provided a new technique for providing attenuation correction in 
emission computed tomography. 
[Solving Means] 

A method of ML-EM image reconstruction is provided for use in connection with a 
diagnostic imaging apparatus (10) that generates projection data. The method includes 
collecting projection data, including measured emission projection data. An initial emission 
map and attenuation map are assumed. The emission map and the attenuation map are 
iteratxvely updated. With each iteration, the emission map is recalculated by taking a 
previous emission map and adjusting it based upon: (5) the measured emission projection 
data; 00 * ^projection of the previous emission msp which is carded out with a multi- 
dimensional projection model; and, (Hi) a projection of the attenuation map. As well, 
with each iteration, the attenuation map is recalculated by taking a previous attenuation map 
and adjusting h based upon: (0 the measured emission projection data; and, (H) a 
rcprojection of the previous emission map which is earned out with the muM-dimenaional 
projection model. In a preferred embodiment, with source-assisted reconstruction, the 
recalc u la t ion of die attenuation map is additionally based upon: (Hi) measured transmission 
projection data; and (pr) a reference or blank data set of measured transmission projection 
data taken without the subject present in the imaging apparatus (10). 



[Selected Drawing] FIG. 7 



